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The drag prediction workshop series (DPW), held over the last six years, and sponsored by 
the AIAA Applied Aerodynamics Committee, has been extremely useful in providing an 
assessment of the state-of-the-art in computationally based aerodynamic drag prediction. 

An emerging consensus from the three workshop series has been the identification of spa- 
tial discretization errors as a dominant error source in absolute as well as incremental drag 
prediction. This paper provides an overview of the collective experience from the workshop 
series regarding the effect of grid-related issues on overall drag prediction accuracy. Exam- 
ples based on workshop results are used to illustrate the effect of grid resolution and grid 
quality on drag prediction, and grid convergence behavior is examined in detail. For fully 
attached flows, various accurate and successful workshop results are demonstrated, while 
anomalous behavior is identified for a number of cases involving substantial regions of sepa- 
rated flow. Based on collective workshop experiences, recommendations for improvements 
in mesh generation technology which have the potential to impact the state-of-the-art of 
aerodynamic drag prediction are given. 

Over the last six years, the AIAA Applied Aerodynamics Committee has sponsored three Drag Prediction 
Workshops (DPW), with the aim of assessing the state-of-the-art of current Computational Fluid Dynamics 
(CFD) solvers at predicting absolute and incremental drag changes on generic transonic transport aircraft 
configurations. The first workshop, DPW I, held in Anaheim CA, in June 2001, employed the DLR-F4 wing- 
body geometry 1,2 due to its simplicity and relevance to modern transonic transport aircraft design studies, as 
well as due to the availability of publicly released experimental data. The second workshop, DPW II, 3 held 
in Orlando FL, in June 2003, employed the slightly different DLR-F6 geometry, which included a wing-body 
configuration, and a wing-body-nacelle-pylon configuration, both of which have been tested experimentally. 
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DPW II thus offered the possibility to study increments in drag between two similar configurations, as would 
be typically done in a design study. Results from both the first and second workshops showed considerable 
scatter in terms of drag prediction and agreement with experimental data, which was largely attributed to 
the fact that the test cases contained substantial regions of turbulent separated flow, which are known to be 
very sensitive to turbulence model fidelity, numerical dissipation, and discretization errors. 

For the third and most recent workshop, DPW III, 4 held in San Francisco CA, in June 2006, a wing-body 
fairing was designed 5 and added to the DLR-F6 wing body configuration in order to suppress flow separation 
in the wing-body juncture region. The baseline DLR-F6 wing-body configuration was also retained as a test 
case, providing a bridge between the two most recent workshops, and enabling the comparison of incremental 
values between clean and separated flow cases. Furthermore, in order to reduce the extent of the trailing 
edge separation observed on the baseline DLR-F6 configuration studied in DPW II, the Reynolds number 
for these cases was raised from 3 million to 5 million in DPW III. In addition to the two DLR-F6 wing-body 
configurations, a set of two wing-alone configurations were included as a second test case for DPW III. The 
purpose of this second test case was to provide a simpler wing-alone geometry, which would enable the use 
of higher grid resolution, thus enabling a more complete grid resolution study, and to provide a well behaved 
test case with no anticipated regions of separated flow. While one of the objectives of DPW III was to 
test the hypothesis that the large scatter in the results of the previous workshops was due to the presence 
of separated flow, considerable scatter was also observed for the DPW III results, even for the cases with 
negligible flow separation such as the wing-body-fairing cases and the simple wing cases. Furthermore, in 
many instances this scatter was not observed to decrease with increasing grid resolution. 

Over the course of these three workshops 6 8 and numerous follow-on studies undertaken over these 
years, 9-27 a consensus has been emerging suggesting that the dominant sources of error in aerodynamic 
drag prediction are related to spatial discretization error, which in turn stems largely from grid resolution 
and grid quality issues. The gradual evolution of this consensus can be seen in the focus of the three DPW 
workshops. In DPW I, the workshop results were all presented for a “standard” grid size, which at the 
time was in the neighborhood of 3 million points. 6 The considerable scatter of the collective results of the 
workshop, along with subsequent independent grid resolution studies, 10 13 prompted the organizers to man- 
date a grid resolution study as part of DPW II. Although the results of the second workshop were analyzed 
more rigorously with regards to spatial error or grid convergence, various anomalies were identified such as 
inconsistent grid convergence behavior, and collective scatter which remained independent of the level of 
grid resolution. The fact that many of these cases contained substantial regions of separated flow further 
clouded the issue of which error source remained most important (modeling or discretization), an issue which 
the third workshop attempted to address through both the addition of the fairing, the inclusion of simpler 
wing alone geometries, and the use of a more extensive grid refinement exercise with much finer grids than 
previously employed. The fact that the bulk of the DPW III results have been presented in the context of a 
grid convergence study is a testament to the importance attributed to determining the effect of discretization 
error on solution accuracy. 

The identification of spatial discretization effects as a dominant error source has important implications 
for advancing the state-of-the-art of aerodynamic simulations, since improvements in other areas such as 
turbulence and transition modeling will either be inconsequential or impossible to assess in the presence 
of dominant spatial discretization error. Thus, understanding and quantifying grid quality and resolution 
effects on solution accuracy remains one of the most important requirements for improving aerodynamic pre- 
dictive ability. The collective DPW results have established some surprising facts concerning the sensitivity 
of aerodynamic drag prediction with respect to computational meshes. Over the course of the workshop 
series, the level of mesh resolution employed for calculations has steadily increased, from grids of several 
million points in DPW I, 6 up to grids of over 40 million points in DPW III, 8 with some individual studies 
employing grids up to 100 million points. 20,22,25 This steady increase in available grid resolution has been 
enabled largely due to advances in computational hardware and improved solution algorithms. However, the 
collective workshop results have not demonstrated substantially reduced variability over the course of the 
workshop series, in spite of over an order of magnitude increase in grid resolution from the first to the third 
workshop. To be fair, much of the observed variability in the collective results stems from the large numbers 
of different solvers, turbulence models, grid types, and grid instances used throughout the workshop series 
by the participants. Individual participant workshop results have demonstrated monotone grid convergence 
behavior, particularly for attached flow cases. However, in many instances, the demonstration of apparent 
grid convergence is contradicted by results obtained with the same code on grids of different construction. 
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These results raise some important questions with respect to the effect of the computational mesh on aero- 
dynamic predictive ability, such as: how can one minimize grid induced simulation errors, what constitutes 
a “good” mesh for aerodynamic simulations, and how should one go about assessing grid convergence ? 

We attempt to address some of these questions in this paper. In the following section, we begin with 
a description of the DPW geometries and test cases, followed by an overview of the gridding guidelines 
established and used for the provided workshop grids. In the subsequent section, we discuss some of the 
important issues related to mesh generation for aerodynamic simulations, such as the establishment of 
adequate grid resolution, the ability to assess grid convergence, and the determination of grid quality on 
solution accuracy. Finally, we suggest a path forward for the community for addressing some of these open 
questions. 


I. Geometry and Test Case Description 

The baseline wing-body configuration used for both DPW II and DPW III, referred to as the DLR-F6, 
corresponds to a transport aircraft wind-tunnel model designed to cruise at transonic speeds. The basic 
planform characteristics of this model are depicted in Figure 1, including the position of the wing mounted 
nacelles, which were part of the DPW II study, but omitted for DPW III. The wing has a leading-edge 
sweep of 27.1°, a quarter chord sweep of 25°, and a dihedral angle of 4.787°. The DLR-F6 model has an 
acute angle of roughly 60° between the fuselage and wing upper-surface at the trailing-edge, which leads to 
flow separation at the wing-body juncture. A wing-body fairing was developed by Vassberg et al. 5 with the 
goal of completely eliminating the separation bubble in this region at the design conditions. The DLR-F6 
wing-body with additional fairing used in the DPW III study is referred to as the DLR-F6 with FX2B or 
simply FX2B configuration. Figure 2 illustrates the computed flow on the DLR-F6 and FX2B configurations 
at the design condition (Mach=0.75, Cl— 0.5, Re=5 million), illustrating the separation bubble for the 
original configuration, and fully attached flow for the new FX2B configuration. DPW III was designed as a 
blind test since no experimental data on the FX2B configuration was available at the time of the workshop, 
although wind-tunnel data has recently been obtained on both configurations in the NASA Langley National 
Transonic Facility (NTF), and will be used for comparison with numerical predictions in subsequent studies. 
For DPW III, participants were required to provide data for both the baseline DLR-F6 geometry and the 
DLR-F6 with FX2B fairing configuration at the following flow conditions: 

• Fixed Cl Single Point Grid Sensitivity Study on Three Grids 

- Mach = 0.75, Cl = 0.5, Re = 5 million 

• Drag Polar on Medium Grid 

- Mach = 0.75, Re = 5 million, Incidence = [-3.0°, -2.0°, -1.0°, -0.5°, 0.0°, 0.5°, 1.0°, 2.5°] 

As mentioned previously, DPW III also included a simpler set of wing-alone geometries for enabling more 
extensive grid refinement studies in the presence of simpler geometries with fully attached flow. Two closely 
related wing configurations were designed by DPW committee members, labeled Wing 1 and Wing 2 (or 
W1 and W2). Both have a common simple trapezoidal planform shape with no leading-edge or trailing- 
edge breaks, and incorporate a leading-edge sweep of 17.2° and quarter-chord sweep of 15.0°. The airfoil 
selected for Wing 1 is representative of a supercritical section found on most modern transonic aircraft, 
producing a nearly flat roof-top pressure distribution with a shock located at approximately 57% chord at 
the design conditions of Mach = 0.76 and Cl= 0.50. The second configuration, Wing 2, was created by 
optimizing the Wing 1 configuration for drag reduction at the design point using only camber and twist 
changes. Additional information regarding the design and definition of these wings is provided by Zickhur 
and Tinoco 3 4 in the workshop presentations on geometry review. For DPW III, participants were required 
to compute the following cases: 

• Fixed Incidence Single Point Grid Sensitivity Study on Four Grids 

— Mach = 0.76, Incidence = 0.5°, Re = 5 million 

• Drag Polar on Medium Grid 

- Mach = 0.76, Re = 5 million, Incidence = [-1.0°, 0.0°, 0.5°, 1.0°, 1.5°, 2.0°, 2.5°3.0°] 
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Figure 1. DLR-F6 wing-body planform definition including standard nacelle configuration used in DPW II. 
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Figure 2. Illustration of separated flow in wing-root junction region for (a) DLR-F6 configuration and (b) 
suppression of flow separation in DLR-F6 with FX2B fairing configuration. 
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These wing- alone configurations were devised purely for computational studies, and no experimental data 
exists for these cases. 

Although the workshop results called for both fixed point grid refinement studies as well as the calculation 
of drag polars on a fixed grid, we will concentrate on the results of the fixed-point grid refinement studies 
for the purposes of this paper, with the intent of focusing on grid-related errors for drag prediction. Note 
that the grid refinement studies are specified at fixed Cl for the DLR-F6 and FX2B configurations, while 
these are specified at fixed incidence for the wing-alone cases. 

II. DPW Gridding Guidelines 

The wide range of scales present in aerodynamic problems requires that even the simplest simulations 
employ highly variable mesh resolution, tailored to the problem at hand. In the absence of fully automated 
adaptive meshing procedures, this means that mesh generation for computational aerodynamics will remain 
somewhat of an art based on practitioner experience and intuition, at least for the near term. In an attempt 
to ensure some level of uniformity over the various grids used for the DPW calculations, the committee 
developed a set of gridding guidelines, which codifies much of the collective experience of the committee 
members in aerodynamic grid generation practices. These guidelines have remained relatively unchanged 
over the course of the three workshops, and thus only those published for DPW III are discussed herein. 
For the latter workshops, the additional focus on grid convergence issues also resulted in the issuance of 
guidelines for generating consistent families of uniformly refined meshes for grid convergence studies. The 
DPW III gridding guidelines are summarized below where C,M,F,X refer to Coarse, Medium, Fine, and 
Extra Fine mesh levels, respectively. 

Gridding Guidelines 

• Boundary Layer Region 

-¥+<[1,11,^] [ C,M,F,X ] 

— Ai ~ 0.0006mm, (Approximate dimensional spacing for Y + = 1) 

— A 2 = Ai, (Two cell layers of constant spacing at viscous wall) 

— Growth Rates 1.25 (Preferably 1.20) 

• Farfield: ~ 100 C re f -lengths away from geometry 

• Local Spacings on Medium Grid 

— Chord wise: 0.1% local chord at Wing Leading Edge and Trailing Edge 
— Spanwise: 0.1% semispan at root and tip 
— Cell size on Fuselage nose and tail: 2% C re f 

• Cells across Wing Trailing-Edge Base: [8,12,16,24] [C,M,F,X] 

• Grid Family 

— Medium Mesh Representative of Current Engineering Drag Predictions 
— Maintain a Parametric Family of Uniformly- Refined Grids in Sequence 

— Grid Size to Grow ~ 3X for Each Level Refinement [ Structured: 1.5X in each I,J,K Direction ] 

— Give Consideration to Multigridable Dimensions on Structured Meshes 
- Sample Size for DLR-F6 Wing/Body: [ 2.7M, 8M, 24M ] [C,M,F ] 


The guidelines specify an outer boundary location of approximately 100 reference chords out from the body, 
which has been found to minimize the effect of the outer boundary condition treatment on computed lift 
and drag values. Local grid spacings of 0.1% chord at critical regions such as leading edge, trailing edge, and 
wing tip and root areas are also specified for the medium level grid, with the appropriate scaling of these 
spacings for the coarser and finer grids. For blunt wing trailing-edge geometries, between 8 to 24 cells across 
the trailing-edge thickness are specified, as experience has shown the high level of sensitivity of the results 
to resolution along blunt trailing edges. The standard practice of resolving the boundary layer down to a 
minimum of = 1 is also specified for the coarsest grid, with appropriately smaller spacings for the finer 
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grids. Additionally, the first two layers of cells at the wall should be of similar spacings, and the growth rate 
of the normal spacing through the boundary layer should not exceed 1.25. For grid convergence studies, the 
specifications call for sequences of grids with uniform refinement in all areas of the domain, resulting in a 
growth of the number of grid points by a factor of 3 for each consecutively finer mesh, with approximate 
grid sizes of 2.7M, 8M, and 24M points for the coarse, medium and fine meshes, respectively. 


Table 1. DPW III Case I (DLR-F6/FX2B) Submitted Grids - Number of Grid Points 



DLR-F6 

FX2B 

Family 

Type 

Fine 

Medium 

Coarse 

Fine 

Medium 

Coarse 

Tinoco 

MB 

27,185,664 

8,080,896 

2,298,880 

27,185,664 

8,080,896 

2,298,880 

SAUNA 

MB 

9,761,201 

4,731,073 

2,551,989 

9,761,201 

4,731,073 

2,551,989 

Extruded 

MB 

28,367,120 

9,343,009 

2,996,626 

28,367,120 

9,329,185 

3,028,420 

Gridgen 

MB 

27,982,776 

8,927,196 

2,739,621 

27,982,776 

9,138,772 

2,842,878 

Sclafani 

OS 

26,892,352 

7,985,236 

2,387,918 

26,969,192 

8,020,348 

2,395,170 

DLR 

HY 

8,535,263 

5,102,446 

2,464,385 

10,305,876 

6,111,664 

2,873,102 

ANSYS 

HY 

18,120,772 

8,038,922 

3,059,189 

20,472,520 

8,272,308 

3,163,605 

LaRC 

UH 

40,014,934 

14,298,135 

5,354,214 

41,069,036 

14,598,610 

5,618,073 

AFLR 

UN 

11,374,451 

3,792,485 

1,492,082 

11,849,212 

3,178,559 

1,640,590 

Embraer 

UN 

24,030,000 

8,320,000 

3,550,000 

24,030,000 

8,320,000 

3,550,000 

STAR 

UN 

- 

12,377,058 

- 

21,509,137 

12,469,599 

8,421,799 

USM3D 

UN 

- 

- 

- 

- 

- 

- 

TAS 

UN 

17,535,215 

9,431,154 

5,399,929 

17,219,535 

9,481,477 

5,422,128 


Table 2. DPW III Case 2 (Wing-l/Wing-2) Submitted Grids - Number of Grid Points 



DPW-W1 

DPW-W2 

Family 

Extra- Fine 

Fine 

Medium 

Coarse 

Extra-Fine 

Fine 

Medium 

Coarse 

Tinoco 

14,811,489 

8,620,123 

4,204,203 

1,602,651 

14,811,489 

8,620,123 

4,204,203 

1,602,651 

Sclafani 

55,014,321 

16,265,909 

4,856,149 

1,442,285 

55,014,321 

16,265,909 

4,856,149 

1,442,285 

DLR 

17,053,510 

10,150,588 

5,288,507 

2,174,364 

16,631,805 

9,910,645 

5,030,379 

1,928,405 

LaRC 

36,956,019 

11,492,625 

4,495,117 

1,818,508 

38,462,630 

11,903,329 

4,658,853 

1,882,672 

Raytheon 

12,748,678 

6,138,245 

2,417,082 

983,633 

12,419,567 

5,963,713 

2,325,884 

947,409 

JFLO 

- 

- 

- 

- 

- 

- 

- 

- 


Tables 1 and 2 list the attributes of the various grids submitted and employed for computations for 
DPW III, while Figures 3 through 7 illustrate the topology and resolution for several of the different baseline 
grids used for the workshop. These tables and figures illustrate the wide variety of grid types used for the 
workshop, which include overset structured meshes, point-matched multi-block structured grids, unstruc- 
tured tetrahedral meshes, and hybrid prismatic-tetrahedral meshes. Only the general characteristics of these 
grids are shown here, and further details of the individual grid submissions can be found in companion 
papers. 22-25,27 

Although all submitted grids were generated according to the DPW published guidelines, specific com- 
pliance with the guidelines of individual submissions could not be monitored, and it remains unclear how 
closely all these submissions adhere to the guidelines. Furthermore, specific characteristics of various grid 
generation packages may result in substantially different grid characteristics based on the same guidelines. 
For example, for overset structured meshes, the topology of the grid system and location and extent of the 
overlap regions must be considered carefully for accurate numerical drag prediction, as discussed by Vassberg 
et al. 28 For unstructured meshes, much higher node counts will result for approaches which do not allow 
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or employ spanwise stretching, especially considering the high level of trailing-edge resolution specified by 
the guidelines. In some cases, the trailing-edge resolution must be relaxed in order to maintain practical 
grid sizes in the absence of spanwise stretching. 23 Additionally, for unstructured tetrahedral mesh solvers, 
the specification of mesh resolution for cell-centered discretizations generally does not correspond to the 
specification for vertex-based discretizations, since the latter approach contains 5 to 6 times fewer degrees of 
freedom on the same mesh as the cell-centered approach. However, because vertex-based discretizations entail 
more complex stencils, establishing the equivalent grid resolution for similar accuracy for both approaches 
has remained an open question. Experience from previous DPW exercises has suggested that matching the 
number of triangles on the boundary surface for cell-centered approaches to the number of vertices on the 
surface for vertex-based discretizations (i.e. a 2:1 resolution ratio on the surface), results in vertex-based 
meshes containing three times more points than the equivalent cell-centered mesh, and produces equivalent 
accuracy for both schemes. However, these findings are far from conclusive, and the situation becomes even 
more complex for hybrid meshes containing mixtures of tetrahedra, prisms, pyramids and hexahedra. 

In summary, the wide variety of grid types and technologies used for the workshop series implies that the 
characterization of these grids by simple node counts and guidelines is approximate at best, and substantially 
different characteristics can be expected between supposedly equivalent grids. On the other hand, one of the 
outcomes of the workshop series is a large data-base of publicly available grids and published computational 
results, which can be used to assess grid quality and perform further numerical experiments. 

III. Overview of DPW Grid Related Results 

Figure 8 depicts the collective results from DPW I, while Figures 9 and 10 show the collective results 
from DPW II and DPW III respectively. While the DPW I results were only examined at a fixed mesh 
resolution, the DPW II and DPW III results were examined as a function of grid resolution, using a grid 
index factor in DPW III defined as the number of grid points raised to the —2/3 power, as a measure 
of average cell size squared, which should produce straight line plots for grid convergence studies based 
on second-order accurate discretizations. Furthermore, in DPW III, the idealized drag coefficient (defined 
as C Dideal = Go ~ Cl 2 /(PIAR), where AR refers to the wing planform aspect ratio) is used instead of 
the actual drag coefficient. By subtracting out the estimated induced drag, the use of the idealized drag 
coefficient represents an attempt to minimize the effect of variations in lift of the submitted results. The 
shear amount of information present in these plots is a testament to the vast numbers of different solvers, 
turbulence models, grid types, and grid resolutions used by workshop participants. Since it is well known that 
different solvers can be more or less sensitive or tolerant to different mesh characteristics, the ideal situation 
for assessing grid suitability would be to run a single grid with all solvers, or conversely all grids with a 
single solver. Unfortunately, the large number of different entries in the DPW data-sets, combined with 
the fact that these workshops represent good-will efforts on the part of the community at large, obviate the 
possibility of performing such types of controlled experiments. Nevertheless, various important observations 
can be made from the DPW workshop results regarding the influence of mesh generation technology on the 
accuracy of aerodynamic simulations. In the following, we discuss separately the issues of grid resolution, 
grid convergence, and grid quality. 

A. Grid Resolution 

It is generally evident from the workshop results that significantly increased grid resolution will be necessary 
for meeting the accuracy requirements of aerodynamic drag prediction. To some degree, a natural progression 
in what is considered practical grid resolution occurs as hardware and solver technologies advance. This is 
reflected in the grid sizes used in the three-successive DPW workshops, which have increased by over an 
order of magnitude in this 5 year period, going from approximately 3 million points in 2001 (DPW I) to over 
40M points in 2006 (DPW III), with some interim and follow-on studies using grids in the vicinity of 100 
million points. 20,22,25 However, given the levels of discrepancies noted in even the latest workshop results, 
accompanied by the grid convergence trends, it is evident that the ability to employ much finer grids, perhaps 
up to 10 9 grid points on a regular basis, would greatly improve the accuracy and confidence in production 
level aerodynamic drag prediction. 

These high grid resolution requirements stem not only from the high accuracy requirement of drag 
prediction, but also from the large range of scales which must be resolved in aerodynamic simulations. For 
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(a) Surface mesh topology 


(b) Field view of mesh in wing region 




(c) Wing tip mesh topology and resolution 


(d) Trailing edge mesh topology and resolution 


Figure 3. Illustration of overset mesh topology and resolution details for DLR-F6 configuration reproduced 
from Sclafani et al. 22 
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(a) General view of mesh 



(b) Leading-edge WB junction (c) Trailing-edge WB junction (d) Cross section at WB junction 



(e) Wing leading edge 


(f) Airfoil sectional cut 


(g) Wing trailing edge 


Figure 4. Illustration of block-structured mesh topology and resolution of Boeing generated grids from refer- 
ence 25 for DLR-F6 configuration 
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(a) Block topology 


111 

(b) General view of mesh 



(c) Leading-edge WB junction 


(d) Trailing-edge WB junction 


(e) Cross section at WB junction 



(f) Wing leading edge 


(g) Airfoil sectional cut 


(h) Wing trailing edge 


Figure 5. Illustration of block-structured mesh topology and resolution of JAXA generated grids from refer- 
ence 23 for DLR-F6 configuration 
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Figure 6. Wing-body-fairing unstructured grid resolution near wing-fuselage juncture for DPW III meshes 
generated by NASA Langley and designated LaRC node-centered grids: (a) coarse mesh (5.6 million points 
total), (b) medium mesh (14.6 million points total), (c) fine mesh (41 million points total). 



Figure 7. DPW III Wing-1 fine unstructured mesh generated by NASA Langley using VGRID software, (a): 
Fine Grid on Wingl Configuration comprising a total of 15 million points, (b): Detail of trailing-edge base 
resolution for fine grid on Wingl Configuration containing pproximately 16 cells across blunt trailing-edge 
base. 
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Multiblock Structured Unstructured 




Figure 8. Collective results from DPW I plotted as computed idealized drag coefficient versus lift coefficient 
for various solvers and grid types. 
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Figure 9. Collective results from DPW II plotted as computed drag coefficient at fixed lift condition versus 
number of grid points (a) wing-body configuration and (b) wing-body-nacelle-pylon configuration. 
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Figure 10. Collective results from DPW III plotted as computed idealized drag coefficient at fixed lift condition 
versus grid resolution index for DLR-F6 wing-body and wing-body with FX2B fairing. 


example, in the vicinity of the outer boundary, the dominant scales are many times larger than the mean 
chord or characteristic body dimension, while flow details in the vicinity of blunt trailing edges, where the 
characteristic length is thousands of times smaller than the mean chord, can have a profound effect on lift 
and drag accuracy if not properly resolved. Additionally, the thin boundary layers produced at test and 
flight Reynolds numbers can only be resolved in a tractable way using highly stretched anisotropic meshes 
with aspect ratios in excess of 1000:1 or even 10000:1. These characteristics dictate the use of highly variable 
resolution and anisotropic meshes for aerodynamic calculations, and are the principal reason mesh generation 
has remained an art largely based on practitioner experience and intuition. 

The highly variable and tailored resolution required for accurate drag prediction is illustrated in Figures 
3 through 7 for representative sets of overset structured, point-matched block-structured, and unstructured 
meshes used in the workshop. For the structured mesh sets, the figures demonstrate how the blocking or 
overset topology can have profound effects on the resulting local resolution. For the unstructured meshes, 
issues such as element type and spanwise stretching can play a significant role in the resulting resolution as 
well. In all cases, the figures illustrate highly variable resolution resulting from attempts to provide enhanced 
refinement in critical areas, such as at blunt trailing edges, and at the wing- fuselage junction, where a pocket 
of separated flow is anticipated for the DLR-F6 configuration. 

The workshop results have demonstrated the complex nature of the interaction between grid resolution 
and drag prediction accuracy. In most instances, increased grid resolution is found to have only a minor effect 
on local quantities such as surface pressure coefficients or skin friction. This behavior is illustrated in Figure 
11, for the W1 test case, using four successively finer unstructured meshes from 1.8 million to 37 million 
points, using the NSU3D unstructured mesh solver. 27 In this case, the only notable differences in the surface 
values in going from coarse to fine meshes are in the resolution of the shock wave, which becomes sharper, 
while rotating about its center point. While in this case the shock position remains relatively unchanged, 
cases where the shock position changes slightly with grid resolution result in shocks of different strength, 
thus altering the computed drag values appreciably. In such cases, changes in the shock position are due 
most often to inadequate resolution ahead of the shock, which thickens the boundary layer through the 
convection of upstream generated numerical entropy. Similarly, in addition to their direct effect on friction 
drag, inaccuracies in computed skin friction can alter the prediction of downstream separation locations, 
which results in a much more profound effect on overall drag prediction. An example of this behavior was 
originally given Anderson et ah, 29 and reproduced in reference, 30 where the effects of inadequate boundary- 
layer resolution were found to produce premature separation on a two-dimensional high-lift configuration. 
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(a) (b) 

Figure 11. (a) Computed surface pressure coefficients and (b) skin- friction coefficients at the 55.1% span 

location at Mach = 0.76, a = 0.5°, Reynolds = 5 million, on the four grids used in the grid refinement study 
for Wingl with the NSU3D unstructured mesh solver (reproduced from Mavriplis 27 ). 

Examination of individual workshop results has shown that grid resolution requirements for drag predic- 
tion involving regions of separated flow are much more demanding and less well understood than for cases 
where the flow is fully attached. For example, various participants have reported adequate grid convergence 
of drag values as well as relative insensitivity of these values to grid topology and grid type for the FX2B 
and wing-alone configurations in DPW III. On the other hand, for the DLR-F6 wing-body case involving 
pockets of separated flow, many participants have observed continuing sensitivity of the computed drag 
values with increasing grid resolution, as well as with differing grid topologies and grid types. Furthermore, 
no consensus on the sensitivity of drag prediction to local variations in grid resolution can be determined 
from the workshop results. For example, Sclafani et al. 22 observe that the size of the separation bubble at 
the wing root junction of the DLR-F6 (and thus the drag value) increases continuously with increasing grid 
resolution for an overset mesh approach, while Tinoco et al. 25 and Murayama et al. 23 observe a relatively 
constant separation bubble size with increasing grid resolution for a block-structured mesh approach. Using 
an unstructured mesh approach, Murayama et al. 23 also observe little variation in the amount of flow sepa- 
ration with increasing grid resolution, while Tinoco et al. 25 observe large increases in flow separation using 
two different unstructured mesh solvers. Similar inconsistencies are observed with respect to trailing-edge 
separation, which is found to increase as a function of grid resolution by Mavriplis 27 using an unstructured 
mesh approach, but found to decrease or remain relatively constant by Murayama et al. 23 for both their 
unstructured and block-structured approaches, and by Sclafani et al. 22 using OVERFLOW. 

In the absence of a clear understanding of the necessary grid resolution requirements, experimentation 
with progressively finer meshes remains the most viable approach. Grid refinement studies are predicated 
on the assumption of asymptotic behavior of the solution error, given some adequate level of grid resolution 
throughout the computational domain. Surveying the DPW III grid convergence study results, one may 
conclude that the coarse baseline meshes used in the workshop, containing of the order of 2 to 2.5 million 
points, often do not lie in the asymptotic region for absolute drag prediction, even for fully attached flow 
cases (see for example references. 8,22 ) Thus, the minimum adequate grid resolution for asymptotic behavior 
appears to be in the range of 5 to 10 million points, although this can also be expected to depend strongly 
on the optimality of the distribution of grid resolution throughout the domain. 

One of the problems with increasing grid resolution is the rapid growth in grid sizes with successive 
refinements. In three dimensions, a uniform grid refinement by a factor of 2 in each direction results in 
an overall increase in grid size by a factor of 8, with at best an equivalent (but often larger) increase in 
required flow solution resources, for an optimally convergent and scalable flow solver. Thus, a single uniform 
refinement of one of the fine DPW III baseline meshes of 40M points would result in a mesh of 320M points, 
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which is already beyond the capabilities of most practitioners. Furthermore, few mesh generation packages 
are capable of generating such large scale meshes in a robust and efficient manner. 

In order to evaluate the suitability of a particular mesh for a given aerodynamic simulation, an a posteriori 
error-estimation technique is required which, when coupled with an adaptive meshing strategy, can be 
used to improve the mesh resolution and accuracy of the simulation. Figure 12 provides an illustration of 
an adapted unstructured mesh used for improving the accuracy of drag prediction in DPW III employed 
by Brodersen et al. 24 In this case, the mesh was refined locally based on the gradient of total pressure 
in the computed solution. The principal difficulty in effective use of adaptive mesh refinement for drag 
prediction problems lies in the reliability of these types of gradient-based error indicators, which themselves 
require the solution to he in the asymptotic error region. More sophisticated approaches such as adjoint- 
based functional error-estimation techniques have shown promise, 31 ~ 33 but these also rely on a linearization 
about the current solution, and will require further development to achieve the level of reliability sought for 
production aerodynamic calculations. 



Figure 12. Illustration of adaptively refined mesh for DPW III Wing-1 drag calculation using gradient of total 
pressure as a refinement criterion, and limiting the total increase in number of grid points to 70%. (a) Surface 
grid (b) Field slice at 77 = 0.17 (reproduced from Brodersen et al. 24 ). 


B. Grid Convergence 

As mentioned previously, the DPW workshop series has gradually devoted increased attention to the issue of 
grid convergence. The principle concept of grid convergence states that as the mesh resolution is increased 
in all areas of the domain, the discretization error in the solution should vanish and the discretized solution 
should gradually tend towards the continuous solution. For a second-order accurate scheme, where the 
discretization error varies as 0(h 2 ), and where h represents some measure of the mesh spacing, solution 
values such as lift and drag plotted versus h 2 should produce a straight-line plot which extrapolates to the 
continuous solution values in the limit as h — > 0 . The DPW III grid convergence studies consist of plotting 
solution values versus a grid index factor, defined as TV - 3 ? where N represents the number of grid points. 
The grid index factor is meant to be representative of an average mesh spacing squared (i.e. h 2 ). However, 
since the grids all contain highly variable mesh spacings, two grids with vastly different local grid spacings 
can have similar grid index factors, and this measure only makes sense when applied to families of self-similar 
grids. 34 A family of self-similar grids consists of a set of coarse and fine meshes where the resolution of the 
finer meshes is obtained by uniformly reducing the mesh spacing of the coarser meshes in all regions of the 
domain. For overset- or block-structured meshes, this can be achieved by specifying increased numbers of 
points in the I,J,K coordinate directions. For unstructured meshes, a global mesh spacing factor is required, 
which is used to rescale the specified spacings in all areas of the domain in a uniform manner. Although 
many unstructured mesh generation packages such as VGRID 35 have been designed with the use of a global 
grid spacing factor in mind, experience has shown that the resulting meshes do not always conform to the 
desired uniform refinement characteristics throughout all regions of the domain, due to specifics of the mesh 
generation algorithm. Additionally, for both structured and unstructured meshes, the algorithmic parameters 
which determine the exponential growth of the normal boundary-layer resolution must be considered carefully 
in order to approximate as closely as possible the goal of uniform refinement of the anisotropic mesh in this 
region. The generation of consistent families of refined meshes for grid convergence studies in the presence 
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of complex geometries is therefore a non-trivial task. The degree to which the uniform refinement objectives 
have been met for the families of grids used in the DPW workshops is not entirely clear, although it can 
be presumed that the efforts of the workshops represent the state-of-the-art in current grid generation 
technology. 

Another issue which arises in the context of a grid convergence study is the ability to generate and employ 
sufficiently fine meshes. For grid convergence plots, a minimum of three data points (and preferably four 
or more points) is required in order to establish the degree of grid convergence of the solver. Since each 
successive uniform mesh refinement by a factor of two results in eight times higher grid resolution, using a 
four-grid sequence requires a factor of 512 increase in the number of grid points between the coarsest and 
finest mesh of the sequence. The situation quickly becomes intractable, particularly given the previously 
stated observation that the coarsest meshes of the sequence must often be larger than several million points in 
order to fall in the asymptotic region. For this reason, the sequence of meshes in the DPW grid convergence 
studies were constructed using a lower refinement ratio resulting in a threefold growth in the number of grid 
points between successively finer meshes (mesh spacing reduction of 1.5:1 versus 2:1). However, it is evident 
from this discussion that the ability to generate and employ much finer meshes remains crucial towards 
improving our ability for establishing grid convergence. 

The specific nature of the test case must also be considered carefully in the context of a grid convergence 
study. The majority of the DPW test cases were computed using a fixed Cl condition, rather than a fixed 
incidence condition, since this is the industry practice for drag prediction. For example, the collective DPW 
III results shown in Figures 10 depict computed (idealized) drag at a fixed Cl = 0.5 as a function of the grid 
index value. However, for a truly grid convergent solution, all values including lift and drag should converge 
to their continuous values as the mesh is refined for fixed geometric conditions (i.e. incidence). For these 
reasons, the simpler W1/W2 studies in DPW III were performed at fixed incidence rather than at fixed Cl 
conditions. 

The DPW results have shown that grid convergence can be achieved by a number of solvers when carefully 
constructed families of uniformly refined grids are used, particularly for attached flow cases. Figure 13(a) 
provides an illustration of the grid convergence achieved for the overset structured grid code OVERFLOW 
on the DLR-F6 and FX2B test cases of DPW III, for the computed drag at constant lift on a sequence 
of grids ranging from 2.4 million to 82.6 million points, reproduced from Sclafani et al. 22 The results for 
the FX2B case clearly display second-order accuracy, as evidenced by the straight-line nature of the plot, 
and can be extrapolated reliably to an infinite resolution value by following this line to its intercept on 
the y-axis. On the other hand, the DLR-F6 wing-body results do not produce a straight-line plot, as the 
drag values initially decrease, and then subsequently increase with increasing grid resolution, due to the 
growing presence of separated flow regions with increased grid resolution, as discussed previously. For the 
DPW III wing-alone test cases, reasonable grid convergence is demonstrated in Figures 13(b) and (c) using 
OVERFLOW, 22 although the drag convergence at fixed lift conditions in Figure 13(c) is seen to be superior 
to that at fixed incidence conditions in Figure 13(b). 

Various other workshop participant results have shown similar trends for all grid types, including overset, 
block-structured, and unstructured meshes. Figure 14 illustrates the grid convergence study for the FX2B 
configuration performed by Murayama et al. 23 using block-structured and unstructured meshes. In practice 
it is strictly not valid to plot both types of grids on the same chart with a single grid index factor, since 
they represent two dissimilar families of grids. However, the results demonstrate that both codes achieve 
close to second-order accuracy on their respective family of grids. Perhaps more importantly, extrapolation 
of these plots also reveals that both codes converge to the same drag values in the limit of infinite grid 
resolution, providing a strong verification of the consistency of these two sets of simulations on different 
grid types. On the other hand, using an alternate discretization scheme (termed U-MUSCL in reference 23 ) 
both codes remain consistent among themselves, but converge to a different drag value. The fact that the 
slope of the line for the unstructured mesh results is steeper than the slope of the line for the structured 
mesh results is indicative of higher errors associated with the unstructured mesh results at a given level 
of resolution. In Figures 14(b) and (c) these drag values are further broken down into their component 
contributions (i.e. drag computed on the wing alone and the fuselage alone) and it is observed that the 
slope of the unstructured grid convergence plots is much lower for the wing alone values than for the body 
alone values, thus suggesting inadequate resolution of the fuselage as a principal source of error in the overall 
unstructured mesh drag values. This component drag decomposition strategy provides a rudimentary but 
effective approach for identifying regions of the computational domain requiring additional resolution. 
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DLR-F6 Wing/Body Results: Total Drag 



(a) 


DPW Wing-Alone Results: Total Drag DLR-F6 Wing/Body Results: Total Drag 




(b) (c) 

Figure 13. Grid convergence results from DPW III using overset structured mesh solver OVERFLOW, (a) 
Drag convergence for DLR-F6 with and without FX2B fairing at constant lift, (b) Drag convergence for 
wing-alone geometries at constant incidence, (c) Drag convergence for wing-alone geometries at constant lift 
(reproduced from Sclafani et al. 22 ). 
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Figure 14. Grid convergence for wing-body fairing configuration using block-structured (UPACS) and unstruc- 
tured (TAS) mesh solvers showing different convergence slopes of both methods, and consistent convergence 
of both methods to similar infinite resolution drag values as indicated by the extrapolation arrows for standard 
and modified (U-MUSCL) discretizations, (a) Drag for complete wing-body configuration, (b) Drag due to 
wing component alone, (c) Drag due to body component alone. 
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The simple wing alone (W1-W2) configurations for DPW III were designed to enable an enhanced grid 
convergence study for cases consisting of mostly attached flow, using very fine meshes, and using fixed 
incidence conditions as opposed to the fixed lift conditions employed in the previous cases. 

Figure 15 (a) displays the grid convergence achieved by the NSU3D unstructured mesh solver 27 on 
this case, for W1 and W2, at fixed incidence, using meshes ranging from 1.8 million to 38 million points, 
demonstrating consistent second-order accurate results, as evidenced by the straight-line behavior of both 
plots. In spite of these encouraging results, it has also been shown that when the same calculations are 
performed on grids of different constructions, conflicting grid convergence results can be obtained. This 
is illustrated in Figure 15(b), where a separate grid convergence study for the DPW III W1-W2 case was 
performed on a different family of uniformly refined grids, leading to a different set of straight-line plots which 
do not coincide with the previous plots. Because these two grid convergence studies make use of different 
families of grids, the use of a single grid index value to characterize both families is inappropriate, and there 
is no reason to expect the plots to overlap using this measure. However, consistent grid convergence behavior 
dictates that these curves should asymptote to the same values in the limit of infinite resolution, which is 
clearly not the case in this experiment as seen from Figure 15(b). 

On the one hand, due to the wide variation in spatial grid resolution, it is entirely possible that some 
regions of the domain remain severely under-resolved, even for the finest meshes of the sequence, and that 
the apparent straight-line grid convergence indicative of asymptotic error behavior will be altered as these 
regions become more highly resolved. In this case, one or both of the convergence plots will eventually 
deviate from the currently observed straight line behavior as the grids are further refined, and results from 
both families will eventually converge together. 

On the other hand, another more troubling possibility is that these two families of grids continue con- 
verging to different values as the grids are further refined, which would indicate that at least one of the grid 
configurations is inconsistent, i.e. produces errors of order 0(1). Obviously, to thoroughly investigate this 
behavior requires the ability to generate and employ much finer grids than currently available. 

Other workshop results for this case, such as those previously discussed in Figures 13(b) and (c), have 
suggested that more consistent grid convergence behavior for drag is obtained when these numerical experi- 
ments are performed at fixed lift conditions. In order to explore this phenomena, the convergence of lift at 
fixed incidence for the DPW III wing-alone cases using the NSU3D unstructured mesh solver is plotted in 
Figure 16(a), where it is observed that the lift values computed on the two different families of grids appear 
to be converging to different values. In Figure 16(b), the idealized drag values computed at fixed incidence 
using NSU3D are plotted for both grid families, and these are seen to converge to the similar values in the 
continuous limit. 

These results indicate that the apparent grid convergence inconsistencies are due principally to the 
prediction of lift values, which in turn affect the drag values through the induced drag. However, as stated 
previously, consistent grid convergence implies that all computed values such as drag and lift ultimately 
should converge at fixed incidence. 

C. Grid Quality 

The above grid convergence studies highlight the importance of considering some measure of grid quality 
for specific aerodynamic simulations objectives. Here we implicitly make a distinction between grid quality 
and grid resolution, noting that a grid of poor quality can be refined extensively and still lead to inaccurate 
solutions. Thus, in a grid convergence study, a grid quality metric may be expected to apply equally to an 
entire sequence of grids from within the same family. 

Although grid quality is not directly related to grid resolution, one aspect of grid quality is related to 
the relative distribution of grid resolution in space. For example, a globally fine grid with limited trailing- 
edge resolution would qualify as a poor quality grid for aerodynamic drag prediction. Other characteristics 
which contribute to grid quality include the grid topology, element type, aspect ratio, and orthogonality 
for structured or semi-structured (boundary-layer) grids. Grid quality metrics may differ substantially for 
different types of meshes. Experience-based guidelines for generating high-quality overset meshes for absolute 
drag prediction have been outlined Vassberg et al. 28 These include careful consideration of grid topology, 
avoiding singular lines in critical regions, as well as locating overlap regions in non-critical areas, and matching 
mesh resolution in overlap regions. For unstructured meshes, solution accuracy may be affected by the choice 
of different element types. Although many unstructured mesh solvers make use of fully tetrahedral elements, 
the most common approach is to employ prismatic elements in the highly stretched boundary- layer regions, 
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Figure 15. Grid convergence plots for drag at fixed incidence for Wingl and Wing2 configurations at 
Mach=0.76, Reynolds = 5 million, a = 0.5° using NSU3D unstructured mesh solver for (a) Original work- 
shop meshes ranging from 1.8 to 37M points, and (b) comparison with second family of grids generated at 
Cessna Aircraft Co. 




(a) 


(b) 


Figure 16. Comparison of both grid convergence at fixed incidence for Wingl and Wing2 configurations at 
Mach=0.76, Reynolds = 5 million, a = 0.5° using NSU3D unstructured mesh solver for (a) Lift coefficients and 
(b) Idealized drag coefficients. 
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which results in semi-structured, nearly orthogonal meshes in these areas. A precise quantification of the 
accuracy effects of using fully tetrahedral meshes versus hybrid tetrahedral-prismatic meshes is still lacking, 
and this effect can be expected to differ for vertex-based and cell-based unstructured mesh discretizations. 
Other types of elements such as unstructured sets of hexahedra may also be employed, with possibly beneficial 
effects for cases with multiple stretching directions, such as cases involving spanwise stretching, since the use 
of hexahedral elements avoids the formation of small angle elements in such situations. Figure 17 depicts an 
unstructured hexahedral mesh employed by Brodersen et al. 24 in the computation of the DPW III DLR-F6 
and FX2B geometries. Comparisons between the results computed on this hexahedral mesh and the more 
traditional mixed tetrahedral/prismatic meshes using the DLR tau code have shown minimal differences in 
predicted drag values for the FX2B test case, although somewhat larger differences were seen for the DLR-F6 
separated flow cases. For structured mesh simulations, different grid topologies in the wing-root junction 
area have been investigated by Murayama et al. 23 The two different mesh topologies are depicted in Figure 
18, with the extruded topology shown in Figure 18(a) corresponding to the traditional approach used by 
most workshop participants. Differences in predicted drag values using these two different structured grid 
topologies where found to be minimal for the attached flow cases of the FX2B configuration, while more 
important differences were observed for the DLR-F6 case with separated flow regions. The insensitivity of 
the attached flow results to grid element type or topology is encouraging, while the differences observed in 
the separated flow cases are not necessarily unexpected, given the increased sensitivity of the results to grid 
effects in general, and the possibly differing resolution in these areas for the different grid topologies and 
element types. 



(a) (b) 


Figure IT. Unstructured hexahedral mesh used in reference; 24 (a) Surface grid, (b) Field cut at 30% span 
location. 



Figure 18. Different grid topologies in wing-root junction area used in reference 23 
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On the other hand, various situations have been documented where different grid constructions have 
lead to substantially different aerodynamic predictions using the same solver, even at high resolution. An 
example of such behavior is illustrated in Figures 19, reproduced from Mavriplis, 20 where the computed lift 
and drag coefficients for the DLR-F6 wing-body configuration from DPW II on a family of uniformly refined 
meshes ranging up to 72 million points are compared with results obtained on an alternate mesh of 65 million 
points. Results on the first set of meshes display consistent grid convergence, going from 1 million up to 
72 million grid points, for both a transonic (Mach=0.75) and subsonic (Mach=0.3) test cases, and values 
computed using two differing values of the dissipation parameter are seen to converge consistently to similar 
infinite resolution values. However, the lift and drag values computed on the 65 million point grid, which 
does not belong to the same family of grids, deviate substantially from the path of grid convergence defined 
by the first set of grids. A 10% drop in lift and corresponding decrease in drag are observed between the 
72 million point grid and the 65 million point grid. Note that it is strictly not valid to plot these two grids 
using the same grid index, since they are not representative of a single family of grids. However, this is done 
in this case to illustrate the abrupt changes between the two finest grids. It remains very surprising and 
unexpected that there could be such a large difference between two solutions computed on such fine grids 
(72M versus 65M pts) especially in light of the apparent grid convergence path determined by the original 
family of grids. It is worth noting that these trends (lower lift and drag) are observed as well for the subsonic 
(Mach=0.3) case, and have been reproduced using the FUN3D unstructured mesh flow solver by Lee-Rausch 
et ah, 36 verifying that these results are not simply the product of anomalous flow solver behavior. Figure 20 
illustrates details of the mesh topology in the blunt wing trailing edge region, where spanwise stretching for 
the 72 million point mesh is evident, while isotropic resolution is observed for the 65 million point mesh. The 
source of the discrepancies in the results computed on these different meshes remains unclear, although one 
may speculate that the use of anisotropic cells in the trailing-edge region results in the prediction of smaller 
trailing-edge separation regions, which in turn produces higher lift values when integrated along the span of 
the wing. On the other hand, the demands of isotropic trailing-edge resolution in the 65 million point mesh 
may result in inadequate resolution for this mesh in other critical regions of the domain, which in turn may 
degrade the overall solution accuracy. 

These results have prompted other follow-on studies using structured and unstructured meshes for this 
configuration. In Sclafani et ah, 22 the effect of additional trailing-edge refinement and lower cell aspect 
ratios (close to unity) was studied using overset structured meshes, but was found to have only minimal 
effect on overall drag prediction, while in Murayama et ah, 23 additional trailing-edge refinement using an 
unstructured mesh methodology was also found to have limited impact on predicted drag values. 




(a) (b) 


Figure 19. Comparison of computed lift (a) and drag (b) coefficients versus the number of grid points to the 
-2/3 power for transonic (M=0.75) and subsonic (M=0.3) cases at 0° incidence on sequence of uniformly refined 
unstructured meshes and alternate mesh of 65 million point for DLR-F6 configuration using the NSU3D solver 
(reproduced from Mavriplis 20 ). 
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Figure 20. Comparison between surface resolution for 72 million point grid (a) and 65 million point grid (b) 
in trailing-edge region in the vicinity of 60% span region illustrating spanwise stretching for the 72 million 
point grid. 


This case highlights the importance and difficulties of developing useful mesh quality metrics. In this 
instance, we have two highly resolved grids which yield substantially differing results, which contrast with 
other grid refinement studies, but which has been verified independently by two unstructured mesh solvers. If 
a correlation between solution accuracy and spanwise cell aspect-ratios could be demonstrated conclusively, 
then a grid quality metric could be established for differentiating between such mesh configurations (although 
it is entirely possible that other grid characteristics are responsible for the observed discrepancies). Many 
metrics have been proposed and used over the years including orthogonality measures (for structured or semi- 
structured grids), dihedral angles and vertex degree for unstructured grids, and min/max ratios of consecutive 
cell volumes, to name a few. However, the relationship between these metrics and overall solution accuracy 
remains tenuous in most cases, and few if any of these metrics have been adopted universally. The situation is 
further complicated by the complex interdependence between flow solver discretization and mesh topology. 
Nevertheless, if particular grid topologies or cell shapes could be shown to produce large or inconsistent 
discretization errors, then anomalous grid convergence results could be predicted and eventually avoided 
through grid quality measures and a priori error estimation. Hence, the development of objective grid quality 
metrics remains an important issue for advancing the state-of-the-art in computational aerodynamics. 

IV. Suggestions for Improvement 

Because discretization error has been identified as one of the leading sources of error for computational 
aerodynamics, improvement in the state-of-the-art will require advances in grid generation technology, in 
addition to the development of new metrics and diagnostic tools. At the same time, the complex interde- 
pendence between grid quality /resolution and flow solver performance and accuracy will dictate increased 
collaboration between the mesh generation and flow solution communities. 

The continued development of robust and efficient high-quality mesh generation strategies remains an 
important objective. However, with increasing resolution demands, capabilities for generating meshes of 
several hundred million mesh points and beyond should be targeted. This may involve increased adoption 
and implementation of parallel mesh generation strategies, as well as more elaborate techniques such as 
in-situ parallel mesh refinement, using coarser template meshes (block-structured or unstructured) in order 
to alleviate disk storage and I/O bottlenecks. 

For grid convergence studies, improved tools for generating consistent families of uniformly refined meshes 
are required, particularly for unstructured meshes. Simultaneously, tools for assessing the compliance of 
generated grids with the necessary grid convergence metrics are required. For example, a grid interpolation 
scheme could be developed where the average local resolution between two similar grids is compared at each 
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cell and where areas which exceed prescribed tolerances are flagged. Such a tool could also be extended to 
evaluate grid compliance with user specified guidelines, such as the guidelines issued by the DPW committee. 
In fact, a good first test case would be to run some of the publicly available grids generated for the DPW 
workshops through such a tool. 

Although prediction of absolute drag remains a difficult problem, accurate prediction of incremental drag, 
such as the change in drag between two closely related configurations, is of the utmost importance for guiding 
design decisions. Although accurate incremental drag prediction is often less demanding than absolute drag 
prediction, the DPW results have shown that for cases with separated flow, even reliable incremental drag 
prediction can be difficult to achieve. In order to obtain high quality increment predictions on two related 
configurations, similar grids topologies and resolution must be maintained as much as possible. Thus, grid 
generation techniques which guarantee minimal changes between configurations are desirable. For cases 
where only shape changes occur, with no changes in geometric topology, mesh deformation techniques may 
be most appropriate for constructing self-similar meshes on configurations of interest. 

The remaining areas of interest can be divided into a priori error estimation, and a posterori error 
estimation. A priori error estimation involves the examination of the mesh prior to flow solution calculations 
in order to estimate the suitability of the mesh for the required calculations. This phase involves the 
development of suitable mesh metrics for predicting solution error. At the simple end, this means computing 
and validating specific metrics, such as aspect ratio, dihedral angle etc. A suite of mesh metrics could be 
built up and published, or kernels for implementation be made available on-line to speed adoption. More 
elaborate a priori error-estimation techniques would include devising numerical tests for specific grids and 
discretizations, such as gradient reconstruction techniques etc. which could be run on the subject grid in a 
single pass and report approximation accuracy. 

Finally, the area of a posteriori error estimation holds great long term promise for improving overall 
solution accuracy, particularly when combined with adaptive meshing strategies. As mentioned previously, 
solution gradient-based error-estimation techniques have met with limited success over the years in this 
area. On the other hand, adjoint-based approaches have shown great promise in this area, and have been 
investigated in simplified settings for some time now. 31 33 In fact, an overview of this approach was given 
at the DPW I meeting in 2001, 37 citing some recent two-dimensional work at the time, 31 and outlining the 
promise this approach holds for aerodynamic drag prediction. However, nearly 7 years later, there is still no 
industrial strength adjoint-based adaptive meshing strategy in production use. This slow rate of progress 
is a testament to the complex issues involved in bringing such ideas to fruition, which include not only 
the mechanics of mesh refinement and dynamic load balancing, but other issues such as adjoint solution 
technology and CAD geometry representation in a dynamic parallel environment. The way forward will 
require increased collaboration between all of these different disciplines. 
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